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csj : 1 Introduction 



As is well-known, there do not exist methods of obtaining the general solutions to most 
of non-linear evolutionary PDEs. Very often in such circumstances the only alterna- 
tive to numerical studies of nonlinear models deliver the symmetry-based methods pQ . 
An important sub-class of the self-similar solutions is formed by the travelling wave 
(TW) solutions. The object of this study is to demonstrate the existence of a number 
of localized self-similar TW solutions, basing on the hyperbolic generalization of the 
convection-reaction-diffusion equation. We pay special attention to the existence of 
solitary waves [2J, shock fronts [3] and some other generalized solutions, such as com- 
pactons and peakons [HOE]. The structure of the study is following. In section 2 
we introduce our model equation and next factorize it to an ODE, describing the set 
of TW solutions. Next we discuss the geometric interpretation of solitary waves and 
outline the way of capturing them. In sections 3 and 4 we perform the local nonlinear 
analysis of the dynamical system, equivalent to the factorized ODE, purposed at stat- 
ing conditions of the wave patterns occurrence. In section 5, we present the results of 
numerical study, revealing the presence of all above mentioned types of TW solutions. 
Finally, in section 6, we discuss the results obtained and outline the ways of further 
investigations. 
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2 Statement of the problem 



We consider the following evolutionary equation (referred to as GBE): 



au tt + u t + uu x - x {u n u x ) x = {u- Ux) ip(u). 



(1) 



Here n, x, U\ are positive constants, a is nonnegative. Equation (CEJ) is a generalization 
of both Burgers equation and the reaction-diffusion equation. Let us note, that the 
term au tt appears when the memory effects are taken into account [8], [9], [TUl E] ■ Some 
particular cases of equation (QQ) were studied in recent years [12j [THl HH EH HH [16] . 
Owing to these studies, the analytical description of a large variety of travelling wave 
(TW) solutions is actually available. 

Present investigations are mainly devoted to the qualitative and numerical study of 
the family of TW solution to GBE. Our aim is to show that under certain conditions 
the set of TW solutions contains solitons, compactons, peakons and some other wave 
patterns. To put it briefly, we maintain the notation traditionally used in more specific 
sense. Thus, soliton is usually associated with the exponentially localized invariant TW 
solution to a completely integrable equations, possessing a number of unusual features 
[2]. Some of these features are also inherited by the compactons [TTJIIH]. We maintain 
the notion to those solutions to (pQ), which manifest similar geometric features as "true" 
wave patterns, known under these names. 

Let us consider the set of TW solutions 



Inserting ansatz (T5]) into the GBE, one can obtain, after some manipulation, the fol- 
lowing dynamical system: 



where A(U) = xU n — aV 2 . By analyzing the factorized system ([3]), we are going to 
formulate the conditions contributing to the appearance of the soliton-like solutions and 
the solutions with compact support, called compactons, and some types of generalized 
TW solutions. Analysis carried out, e.g. in [51 E] shows, that homoclinic trajectories 
bi-asymptotic to saddle points correspond to both soliton-like and compacton-like so- 
lutions. In the first case the homoclinic loop is bi-asymptotic to a simple saddle, hence, 
the "time" which is necessary to penetrate such trajectory is infinite. The closed loop 
representing the compacton is bi-asymptotic to a topological saddle. As a result, the 
"time" of penetration is finite. In fact, the compacton is is a compound generalized 
solution. Its compactly supported nonzero part corresponds to the closed loop, while 
the rest corresponds to the stationary point. 

In order to "capture" the homoclinic trajectory among the other solutions to the 
system ([3]), we are going to state the condition, which guarantee the stable limit cycle 
appearance. Choice of such strategy is based upon the well-known fact that the growth 



u(t, x) = 17(0 = U(x-Vt). 



(2) 



A(U) U = A(U) W, 



A(U) W = (U - V) W - xnU n ~ l W 2 - <p(U) (U - U x ) 
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of the radius of the limit cycle in presence of a nearby saddle point most often leads to 
the homoclinic bifurcation. Application of this prescription to the system ([3]) occurs 
to have some peculiarities, which are characterized below. 

1. The most natural parameter of the bifurcation is the wave pack velocity V. Yet 
its change causes the movement of the line of singular points A (£7) = (singular 
line for brevity) in the horizontal direction. As will be shown below, the presence 
of the topological saddle in most cases is due to the fact that the saddle point 
belongs to the singular line. On account of this, the problem of " capturing" the 
compacton-like solution becomes more complicated, for one must "synchronize" 
the moment of the homoclinic bifurcation with the passage of the singular line 
through the saddle point. 

2. Presence of the singular line delivers an extra mechanism of the limit cycle de- 
struction, competing with the mechanism based upon the homoclinic loop for- 
mation. 

3. When the far end of the limit cycle approaches the singular line just at the 
moment of the homoclinic loop formation, the latter becomes flat and reminds 
triangle. Such loop corresponds to a different type of the generalized solutions, 
called peakons [3 [6] . 

4. The singular line is unmovable when a = 0. In this case, corresponding to the 
parabolic-type model, the analysis becomes much more easy, 



3 Andronov— Hopf bifurcation in the system (Sj) 

The homoclinic bifurcation can occur in the system ([3]) if it has an extra stationary 
point. The function 

<p(U) = (U- U ) m V(17), < U < U 1} 

considered throughout the remaining part of the work, assures the required geometric 
configuration, providing that ip(U) does not change sign within the segment [Uq, Ui]. 

To formulate the conditions which guarantee the limit cycle appearance in vicinity 
of the stationary point (U\, 0), let us consider the Jacobi matrix 



Ji 



A(^) 
-<f(Ui) Ut-V 



In order that (Ui, 0) be a center, the eigenvalues of J\ should be pure imaginary. This 
is so if the conditions 

Trace J 1 = U 1 -V = 0, (3) 

Det J x = A(l7i) (fA - U ) m > (4) 
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are fulfilled. The first condition immediately gives us the critical value of the wave pack 
velocity V cri = U\. The second one is equivalent to the statement that both A(£7"i) and 
ip{U\) are nonzero and have the same signs. 

The next thing we are going to do is a study of the stability of the limit cycle. As is 
well known [201 EI] , this is the real part of the first Floquet index JiCi that determines 
the stability of the periodic trajectory. Depending on the sign of A(£7i), there are two 
possibilities. If A(C/i) > when V = V cri or, in other words, the horizontal coordinate 
of the singular line A(£7») = 0, corresponding to the critical value of the parameter V, 
satisfies the inequality 



then conditions ip(Ui) > 0, and 5fCi < should be fulfilled. In case when [/* (V cri ) > 
Ui, these parameters should have the opposite signs. 

To obtain the expression for 3?Ci, the standard formula contained e.g. in [21 J can 
be directly applied, provided that our system is presented in the following form: 

z\ \ _ ( -Q \ ( zi \ + / F(zi,z 2 ) \ /gs 



z 2 J \£l J \ z 2 J \ G(z 1} z 2 

where Q = JJi ■ u, /x = |A(£/i)|, v = \<p(Ui)\, F(z\, z 2 ) and G(zx, z 2 ) stand for nonlinear 
terms. In this (canonical) representation 3?Ci is expressed as follows [2T] : 

16 3? G\ = Fin + Fi 22 + Gn2 + G 222 + — {Fi 2 (Fn + F 22 ) — 

—G\ 2 (Gn + G 22 ) — F\\ G\\ + F 22 G 22 }. (6) 

By Fijk, Fij we denote the coefficients of the function's F(zi,z 2 ) monomials ZiZjZk, 
zi Zj correspondingly. Similarly, indices Gij Gijk denote the coefficient of the second 
and third order monomials of the function G(z\, z 2 ). 

A passage to the canonical variables (zi, z 2 ) can be attained by the unified trans- 
formation. If the relations (J3j) — (SJ) are satisfied, then, rewriting ([3]) in the coordinates 
yx — U—Ux, y 2 = W, we get the following system: 



where 



+ 1 if A(C/i) > 0, 
-1 if A{Ui) < 0, 



$i(l/i, y 2 ) = KnUr 2 yiV2 (u x + J +0 (|i/ 4 | 4 ) , 

$2(2/1, y 2 ) = -\ v\ {mUi) + Wi) Vi) - *nUr 2 yl {U x + (n - 1)^) + Vl y 2 + O (\y t \ A ) 
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In order to pass to the standard representation, we apply the change of coordinates: 

zi ) = ( V ( m 



This gives us the system (J5J), with 

rri 1 - 2 r r7 - 1 1 

+ 0(N 4 ), (9) 



nnur 2 
F(zx, z 2 ) = z 1 z 2 



y/uUi-e - Zi 



G(z 1 ,z 2 ) = [ € 0(U 1 )z 1 -2ip(U 1 )V^-ey^z 1 z 2 
xnU^ 2 2 



n Z2[f/ iv ^-e(n-l)^i]+0(|^| 4 ). 
The real part of the Floquet index is easily calculated from Qj: 

HCi = Gia (Gn + G 22 ) = -e {\^(Ui)\ v{Ui) + xn|^)| t/r 1 } . 

^ 1 (10) 
The above can be summarized in the form of the following statement. 

Theorem 1. J/A(E/i) • ^(E/i) > 0, and 

|A(Z7i)| <p(Ux) + xnb(t/i)| f/r 1 > 0, (11) 

i/ien m vicinity of the critical value of the wave pack velocity V cri = U\ a stable limit 
cycle appears. 

4 Study of the stationary point ((To, 0) 

Rewriting in the coordinates X = U — U , W we obtain the system 
A(U + X)X = A(U + X) W, 

A(U Q + X)W=[(U 1 -U )-X] X m ^(U + X)- (12) 

-xn (U + X) n ~ x W 2 + (U - V + X) W, 

where A(Uq + X) = x\{Uq + X) n — aD 2 ]. Our aim is to determine the conditions 
ensuring that the stationary point X = W = is a topological saddle, or, at least, 
contains a saddle sector in the right half-plane. The standard theory [22] can be applied 
for this purpose. Our system can be written down in the form 



_±{ X \ = f A(U ) \(X 
dT\W \ A U -V \W 



+ nonl. terms, (13) 



where £ = A(U + X)± 

A = 



(Ut-Ua)^), if m 
0, if m > 2. 
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The linearization matrix of the system (fl3j) is nonsingular if m = 1 and the point U* 
lies outside the segment [U , Ui}. In this case the stationary point (0, 0) is a simple 
saddle and the homoclinic trajectory corresponds to the solitary wave solution. Out 
of this case, the Jacobi matrix has at leas one zero eigenvalue. To study the behavior 
of dynamical system in vicinity of a degenerated stationary point, we use the results 
from [22J. Since we are interested in the case when the trace of the Jacobi matrix is 
nonzero, the analysis prescribed in Chapter IX of [22J is the following. 

1. Find the change of variables (U, X, T) i— > (x, y, r) enabling to write down the 
system ffl3|) in the standard form 

■j- = p 2 ( x , y) > 

GET 

7^ = V + Q2 0, y) , 

GST 

where P 2 y), Q2 y) are polynomials of degree 2 or higher. 

2. Solve the equation y + Q 2 (x, y) = with respect to y, presenting the result in 
the form of the decomposition y = a\ x^ 1 + a 2 x m + .... 

3. Find the asymptotic decomposition 

P 2 (x, y(x))=A m x m + .... 

4. Depending on the values of m and the sign of A m , select the type of the complex 
stationary point, using the theorem 65 from |22j. 

5. Return to the original variables (U, X, T) and analyze whether the geometry of 
the problem allows for the homoclinic bifurcation appearance. 



So let us present the results obtained for the system ffl2|) . First we assume, that m > 
1, the statements of the Andronov-Hopf theorem are fulfilled and the point U* satisfying 
the equation A(U*) = lies outside the segment [U , Ui], when the parameter V reaches 
the second bifurcation value V cr2 , corresponding to the homoclinic bifurcation. Here 
we have three possibilities. 

• U* > U\ when the Andronov-Hopf bifurcation occurs (i.e. V = V cri ). This 
inequality does not change up to the homoclinic bifurcation, when V = V cr2 > 

• The inequalities Uq < U* < U\ take place when V = V cri and the inequality 
changes for U* < Uq when V belongs to a small neighborhood of V cr2 < V cri . 

• [/* < Uq when the Andronov-Hopf bifurcation occurs (i.e. V = V cri ). This 
inequality does not change up to the homoclinic bifurcation, when V = V cr2 < 
V cri . 
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Remark 1. Let us note that in the first case the functions A(U) and ip{U) are 
negative when U G \Uq, U-\\. In the third case both of the functions are positive within 
the given interval. In the second case ip(Uo) is positive, and the factor A(Uo) changes 
the sign from negative to positive as XJ* becomes less than U$. 

For m > 1, the canonical system is obtained by the formal change (X, W) \— > (x, y), 
and passage to the new independent variable r = (Uq — V) T, in each of the above 
cases. As a result of such transformation, we get the following system: 



d Ir = y - {*v - Hn Wo + x ) n " 1 y 2 + xm [(^ - ^o) - A MUo) + f(U ) x + ....]} = 

= y + Q 2 (x, y). 



Presenting y in the form of series y — a\ x^ 1 + a2X^ 2 + ... and solving the equation 
y + Q2(x, y) = 0, we obtain 



Fulfillment of the statements of the Andronov-Hopf theorem implies that A(Z7i) ip{U\) > 
when V = V cri . Previously we assumed that function ip(U) does not change sign 
within the segment [U , Ui}. And this is suffice to conclude that the product remains 
positive, when the parameter V attains the value V cr2 , corresponding to the homo- 
clinic bifurcation. It is quite evident for the cases one and three, because the line 
A([7) = remains on the same side of the segment [Uq, U\]. In the case 2 the situa- 
tion is somewhat different, because A(Uq) is negative for V = V cri , while the function 
ip(Uo) is positive and remains so when the parameter V changes. But the singular line 
A({7) = is located to the left from the point (Uq, 0) when V becomes close to V cr2 , 
and then, in accordance with the Remark 1, the product A(Uq) ip{Uo) is positive, when 
the homoclinic bifurcation occurs. So the coefficient A m in the decomposition (|T6|) is 
always negative. Basing on the classification given in Ch. IX of [22J, it is possible to 
formulate the following statement. 

Proposition 1. Let the statements of the Theorem 1 be fulfilled and the singular 
line A(U) = lies outside the segment [Uq, Ui] of the horizontal axis. Then, for 
m > 2 the origin of the system (T$ ) is a topological saddle, having a pair of outgoing 
separatrices tangent to the vertical axis and the pair of incoming ones tangent to the 
horizontal axis, when m = 2 k, k = 1,2, .... For m = 2 k + 1, k — 1,2,3..., the 
stationary point is a saddle-node with two saddle sectors lying in the right half-plane. 
Two outgoing separatrices of the saddle sector are tangent to the vertical axis while the 
incoming one is tangent to the horizontal axis. 



d x 
dr 



A(l/ +:r) 
U -V 



y,= P 2 (x, y), 



(14) 




(15) 



(16) 
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In the following, we present the analysis of system's (fl2|) behavior in vicinity of 
the stationary point (0, 0), assuming that A(U ) = when V = V cr2 . The results of 
the study occur to depend on whether or not Uq is equal to zero. But both of these 
cases can be analyzed simultaneously. We begin with the case m > 1, for which the 
canonical system is obtained by the formal change (X, W) \— > (x, y) and passage to 
the new independent variable r = {Uq — V) T . As a result, we get the following system: 

dx _ k V^n n\ T T n-k k _ p / \ 

dr ~ u -v l^k=i k\(n-ky. u o x y>— r 2\x, y), 

^ T = y-y^j- o {xy-xn (x + U ) n ~ l y 2 + x m [x - (U t + U )) [i/>{U ) +if/{U )x + ...]} 
= y + Q2(x, y). 

(17) 

Presenting y in the form of series y — a\ x^ 1 + a 2 x^ 2 + ... and solving the equation 
y + Q 2 (x, y) = 0, we convince that the first term of the asymptotic decomposition 
y(x) coincides with ( jT5l) . This is not surprising, since the second equations of the 
systems ( TH1) and (JTTJ) are identical. Inserting the function ( fl5l) into the RHS of the 
first equation of system (FlTI) . we get 



P 2 {x, y(x)) 



-nxf/ "-V(f/o)(gf^x m+1 + ...., if U >0, 
-x^(0)^ x n+m + ...., if U = 0. 



In the case m = 1 a passage to the canonical system is attained by means of the 
transformation 

x = X, y = W + BX, t = (U -V)T, 

where 

B - *™ < 18 > 

In the new variables, our system reads as follows: 



S = u^v (V ~ B v) ELi wh)Ur k * k , = K(x, y), 

f T =y+vkr \ kB2 U o~\ n + !) + B + Huo)] x 2 + .... = y + Q 2 (x,y). 
Solving equation y + Q 2 (x, y) = 0, we obtain 

V = tj^-tt [xB 2 Ur\n + 1)+B + if>(uo)] x 2 + .... 

Uq — V 

Inserting function y(x) into the RHS of the first equation, we finally get 



(19) 



P 2 (x, y(x)) 



- >( U^{Uq)^-^x 2 + ... if U o >0, 



-x^(0)^x n+1 + ... if U = 0. 
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Figure 1: Phase portraits of the system ffl2l) in vicinity of the origin for different values 
of the parameter m, in cases when [/* (V cr2 ) — U 
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So the coefficients of the lowest monomials of the decomposition of P2 (x, y(x)) are 
negative and the following statement, based on the classification given in [22], holds 
true for Uq = 0. 

Proposition 2. 

1. If m > 2, and m + n is an odd natural number, then the origin of the system 
|7?P is a topological saddle having a pair of outgoing separatrices tangent to the 
vertical axis and the pair of incoming separatrices tangent to the horizontal axis. 
For even m + n, the stationary point is a saddle-node with two saddle sectors 
lying in the right half-plane. Its outgoing separatrices are tangent to the vertical 
axis while the incoming one is tangent to the horizontal axis. The nodal sector 
lying in the left half-plane is unstable. 

2. If m = 1, and n is an odd number then the origin of the system corresponding to 
( fiPj) is a topological saddle identical with that of the previous case. For even n, 
the stationary point is a saddle-node identical with that of the previous case. 



Below we formulate the analogous result for Uq > 0. 

Proposition 3. If m = 2k, k = 1,2,..., then, the origin of the system (Tf\) is a 



topological saddle having a pair of outgoing separatrices, tangent to the vertical axis and 
the pair of incoming ones, tangent to the horizontal axis. Form = 2k + l, k = 1,2, 
the stationary point is a saddle-node with two saddle sectors lying in the right half- 
plane. Its outgoing separatrices are tangent to the vertical axis while the incoming one 
is tangent to the horizontal axis. The nodal sector lying in the left half-plane is stable. 
For m = l, and arbitrary n 6 N, the origin of the system corresponding to ( T7Pj) is a 
saddle-node identical with that of the case m > 1. 

The last case we are going to analyze is that with m = ~. The motivation for such 
a choice will be clear later on. In order to be able to apply the analytical theory, we 
rewrite the system ffl2l) . introducing new variable Z = yX: 

f| = f [A(U ) + nU^ 1 Z + ... + Z 2n ~ l ] W, 

% = w [(u -v) + z 2 } - ku [u + zT' 1 w 2 + z [{U x - U ) - Z 2 ] i,(U Q + Z 2 ). 

(20) 

Let us consider the Jacobi matrix of the system (120]) . corresponding to the stationary 
point (0, 0): 

fA(£/ ) 
^UoHUi-Uo) (U -V cr2 ) 



J 



(21) 



Analysis of the matrix ( 12T1) shows, that the origin is a simple saddle when U* (V cr2 ) lies 
outside the segment [Uo, U\}. In the case when U* (V cr2 ) = Uq, one of the eigenvalues 
of the Jacobi matrix is zero, and, in order to classify the stationary point, we use the 
results from [22J. A passage to the canonical variables is attained by means of the 
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transformation 



x = Z, y = W + BZ, t 



(Uo-V)T, 



B 



Ui-Uo 
U -V 



The system resulting from this is as follows: 



ft = W^v) [ n K' 1 X + ... + x 2 - 1 ] {y-Bx) = P 2 (x, y), 

^ = y- w ± m {Bf[nUr 1 x + ... + x 2n - 1 ] (y-Bx)-x 2 (y-Bx)- (22 ) 



Solving the equation y + Q 2 (x, y) with respect to y we obtain a series y = a 2 x 2 + 
.... An outlook of this series' coefficients proves to be unimportant, because they do 
not contribute to the lowest term of the asymptotic decomposition of the function 
P 2 (x, y{x)), which is as follows: 



Let us formulate the result obtained as the following statement. 

Proposition 4.1f m = 1/2, then, in case, when U* (V^. 2 ) lies outside the segment 
[Uq, Ui], the origin of the system |Hj) is a simple saddle. In case when U* (V cr2 ) = 0, 
the origin of the system [2^1 is a saddle-node, having two saddle sectors lying in the 
right half-plane. Its outgoing separatrices are tangent to the vertical axis while the 
incoming separatrice is tangent to the horizontal axis. 

The crucial fact appearing from this analysis is that the stationary points (0, 0) of 
the canonical systems (IT7j) . (IT5|) and (T2"2"j) . depending on the values of the parameters 
m, n, are either saddles or saddle-nodes with the saddle sectors placed at the right 
half-space. The return to the original coordinates does not cause any change in the 
position of the saddle sectors, but it changes the orientation of vector fields □ and the 
angles, at which the outgoing separatrices leave the stationary point. The local phase 
portraits corresponding to the distinct cases are shown in Fig. [TJ reconstructed on the 
basis of the analysis of the relation between (TlTj) - (TT9]) and the system ( TT2l) . 

Returning to coordinates (X, W) in the case when m = 1/2, we obtain the patterns 
of the phase trajectories shown in Fig. [2j The difference between the case when 
U* (V cr2 ) = Uq and when U* (V cr2 ) lies outside the segment [U , Ui] is not essential - 
in both cases the incoming and outgoing separatrices enter the origin tangent to the 
vertical axis. Yet, since in the second case the upper sepratrice is a mirror image 

3 We assume that Uq — V cr2 < 0, otherwise the presumable homoclinic loop will not correspond to 
the wave of " compression" . Fulfillment of this condition is verified during the numerical simulation 
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of the lower one, then one can expect the compactly supported solution to be more 
symmetric. 

Before we start to discuss the results of numerical study of the system ([121) . let us 
analyze to what type of solitary waves will correspond the homoclinic loops, presum- 
ably appearing in system ([3]). In order to analyze this issue, we are looking for the 
asymptotic solution W(X) = ax X Ml + a 2 X M2 + ... of the equation 

A(U + X)W^ = G = 

= [(U - V) + X) W - xn (U + X) n ~ l W 2 - ^(U + X) [(U - U x ) + X] X m , 

(23) 

which is equivalent to the system (Tl2l) . Let us start with the case Uq = 0, for which 
the equation (|23l) can be re-written as follows: 

= X 1+ ^ («! + a 2 X Ml + ...) -VX^ (a x + a 2 X Ml + ...) - (24) 
-xnX n+2 ^~ l (a\ + 2axa 2 X" 2 + ....) - [X m+1 + Ux] X">(0) + .... 

The procedure of solving ( 1241) is pure algebraic: we collect the coefficients of different 
powers of X and equalize them to zero. The lowest power in the RHS is either X Ml or 
X n+2ftl_1 . The number n + 2 fix — 1 cannot be less or equal to fix, because it involves 
the inequality < fix < 1 — n, which is impossible for any natural n. On the other 
hand, if n + 2 /^i — 1 < m, then fix becomes an "orphan" and ax should be nullified. 

So x ra+2Ml_1 cannot be the lowest monomial. From this, it immediately appears 
that the only choice leading to a nontrivial solution is fix = m. 

For Uq > the equation is as follows: 

xX 2 ^ 1 (ax + a 2 X^ 1 + ...) (fix «i + fi 2 a 2 X" 1 + ...) = 

= i>(U )X m [(U 1 -U )-X)- 
-ycn [U^ 1 + ...] (a\ X 2 ^ + ...) + [X + (U - V)] [ax X^ + ...]. (25) 

Using the analogous arguments as before for the present (more simple) case, we con- 
clude that fix = rn. 

The decomposition obtained can be used to asymptotically integrate the equation 

^ = W = axX m + .... 
which immediately gives in the lowest order the expression 

[ ai (l_ m )(£_£ )]T^ ; if m ^l 9 

X 

Cexp[ai(£-£ )], if m=l, 

from which we conclude that the trajectory reaches the origin in "finite time" if m < 1. 
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Unfortunately, in case when m G N + , presented above result cannot be attributed 
to both of the saddle sector separatrices, forming the closed loop. In fact, the incoming 
separatrix of the stationary point (0, 0) in all cases considered here is tangent to the 
vertical axis and therefore cannot be described by the formula ( 1261) when m > 1. The 
above formula, then, describes the asymptotic behavior of the "tail " of the solitary 
wave, corresponding to the homoclinic solution. 

To complete the analysis, we resort to some arguments, concerning the other end of 
the homoclinic trajectory. It is the common feature of almost all the cases considered 
here, that the incoming separatrice is tangent to the vertical axis, with the exception 
of the case when m = 1 and U* {V cr2 ) lies outside the interval (Uq, Ui). This enables 
us to assume that in vicinity of the origin W = —B U a + o (U a ), where < a < 1 and 
B > 0. An approximate equation describing the first coordinate of the separatrix is, 
then, as follows 




= -BIT 7 + o{U u ). 



Performing the asymptotic integration, we obtain, up to the change of notation, the 
solution identical with (I26I) , concluding from this that in all analyzed above cases with 
m, n G iV + , the incoming separatix reaches the origin in finite "time". 

It is obvious, that similar arguments can be applied to the analysis of presumable 
homoclinic trajectory, corresponding to m = 1/2. It appears from the above analysis, 
that the incoming and outgoing separatrces are tangent to the vertical axis and both 
of them reach the origin in finite "time". Let us note that when [/* {V cr2 ) lies outside 
[U , £7 2 ], the asymptotic behavior of the spearatrices is described by the formula 

X = ± ^7 (f-fo), 7 > 0. 

So this is the only case when the homoclinic trajectory corresponds with certain to the 
compactly supported solution of the equation ([!]). 



5 Results of numerical simulation 

Since the parabolic case was discussed in detail in our previous work [19] , we mainly 
concentrate here on the hyperbolic case, corresponding to a > 0. Numerical simulations 
of the system (TT21) were carried out with k = 1, XJ\ = 3, = 1. The remaining 
parameters varied from one case to another. The results of qualitative study evidence 
that, within the variety of parameters that were analyzed, the only case when when 
the homoclinic loop corresponds to the compactly supported solution is that with 
m = 1/2. We discuss the results concerning the details of the phase portraits in terms 
of the reference frame (X, W). The numerical experiments show that in case when 
m = 1/2, n — 1 and U* (V cr2 ) < Uq, both of the separatrices forming closed loop enter 
the stationary point (0, 0) tangent to the vertical axis, Fig. [3l the left picture. The 
right picture shows the compactly supported solution to Eq. Qj, corresponding to the 
closed loop. When U* {V cr2 ) tightly approaches Uq, no matter from the left of from the 
right, the left side of the homoclinic trajectory is clasped to the vertical axis, Fig. m 
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Figure 3: Homoclinic solution of the system (fl2|) with (p(Uo + X) = X 1 / 2 (left) and 
the corresponding compactly supported TW solution to Eq. ([I]) (right), obtained for 
n = 1, a = 0.12, Var 2 = 2.68687 and U* - U = -0.133684 
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Figure 4: Homoclinic solution of the system ( fl2|) with tp(Uo + X) = X 1 / 2 (left) and 
the corresponding TW solution to Eq. p]) (right), obtained for n = 1, a — 0.13827, 
V cr2 = 2.68892 and U* - U = 0.99973 



left. The corresponding compactly supported solution has more sharp front, Fig. [U 
right. 

As it was shown in the previous section, appearance of the soliton-like solutions 
is possible merely in case when m = 1. In accordance with the predictions of our 
qualitative analysis, such solutions were observed for m = n = 1 and values of the 
paremeters guaranteeing the fulfillment of the inequality U* (V cr2 ) < Uq. As it is seen 
on the left picture of Fig. [51 both of the separatrices of the saddle sector form nonzero 
angle with the vertical axis and this gives us the reason to state that the homoclinic 
loop corresponds in this case to the smooth solitary wave, which is nonzero for any 

i e R. 

When U* tightly approaches Uq, the incoming separatrice of the homoclinic tra- 
jectory is clasped to the vertical axis (Fig. [6], left) and the solutions reminding shock 
waves with relaxed tails are observed in place of soliton-like wave packs (Fig. [6], center, 
right). 

For m > 2 merely the shock-like solutions have been observed in numerical experi- 
ments. This confirms the arguments put forward in the previous section. The outlook 
of the TW occurs to depend on the values of the parameters m, n. A series of shock-like 
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Figure 5: Homoclinic solution of the system (fl2|) with ip(Uo + X) = X 1 (left), the 
corresponding tandem of well-localized soliton-like solutions to Eq. ([1]) (center), and 
the soliton-like solution (right), obtained for n = 1, a — 0.06, V cr , 2 = 2.65795 and 
U m - U = -0.576119 



4.0 r 4.0 




Figure 6: Homoclinic solution of the system (fl2"j) with tp(Uo + X) = X 1 (left), the 
corresponding tandem od solitary wave solutions to Eq. (TjQ) (center) and a single 
solitary wave solution (right), obtained for n = 1, a = 0.142, V cr2 = 2.65489 and 
U*-U = 0.000878617 
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Figure 7: Tandems od shock-like solutions to Eq. corresponding to (p(Uo + X) = 
X 1 , U* « U , n = 2 (left), n = 3 (center), and n = 4 (right) 




Figure 8: Shock-like solutions to Eq. ([I]), corresponding to <p(U + X) = X 3 , U* ~ U , 
n = 1 (left), n = 3 (center), and n = 4 (right) 



solutions corresponding to m = 1 and n = 2, 3, 4 are shown in Fig. [71 It is seen, that 
effective width of the TW grows as n grows, and the shape becomes more and more 
gently sloping. 

The next series (Fig. [8]) shows the TW corresponding to m = 3 and n = 1, 2, 3. 
This picture differs from the previous one in that the "tails" of the travelling waves 
are longer. A common feature of all the cases with m > 2 is the strong stability of the 
equilibrium, especially in the direction of the outgoing separatrice. Choosing the initial 
data more and more close to the origin, we are able to elongate the "tail" without limit. 
Besides, the profiles of the TW become more and more smooth as the n grows. 

Fig. [9] shows the series obtained for n = 4 and m varying from 1 to 3. It is seen 
that the waves become more and more localized as the parameter m grows. This series 
mainly characterizes the growing stability of the saddle point (0, 0). 

Finally, let us discuss what happens when U* (V cri ) > U\. The birth of the stable 
limit cycle is possible in this case, if <p(U) < 0. Therefore we use in numerical ex- 




Figure 9: Shock-like solutions to Eq. ([T]), corresponding to (f(U + X) = X m , ~ U , 
n = 4, m = 1 (left), m = 2 (center), and m = 3 (right) 
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Figure 10: Periodic solution of the system (fl2l) with (p(Uo + X) = —X 1 (left) and the 
corresponding tandem of generalized cusp-like solutions to Eq. ([1]) (right), obtained 
for n = 1, a = 0.552, V cr2 = 3.00593 and U*-U x = 1.98765 
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Figure 11: Homoclinic solution of the system (Tl2l) with <p(Uo + X) = -I2 (left) and 
the corresponding tandem of generalized cusp- like solutions to Eq. (pQ) (right), obtained 
for n = 1, a = 0.562, V cr2 = 3.14497 and U*-U x = 2.55863 
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periment the function ip(Uo + X) = —X m . Numerical study shows, that the radius 
of the limit cycle grows as the bifurcation parameter V grows. Simultaneously the 
coordinate U*(V) of the singular line intersection with the horizontal axis moves to the 
right, but not so quickly as the radius of the limit cycle. In most cases the destruction 
of the periodic trajectory observed was due to its interaction with the line A(U) = 0. 
For m = n = 1 we succeeded in observing how the limit cycle approaches the line of 
singularity, attaining the triangle shape Fig. [TOj left. The corresponding succession of 
the cusp-like travelling waves is shown in the right picture. Let us note, that we were 
not able to follow the moment of the homoclinic bifurcation, varying the parameters a 
and V. 

Experimenting with m = 1/2, n = 1 occurs to be more easy. The homoclinic 
bifurcation, shown in the left picture of the FigflQ| takes place when the far end of 
the cycle attains the line A(U) = 0. The presence of the singular line causes a drastic 
change of the shape of homoclinic loop, which equally well can be called the "homoclinic 
triangle" . The right picture shows the corresponding succession of the well-separated 
peakons. 

Let us note in conclusion that numerical experiments performed for the parabolic 
case, i.e. for a = 0, demonstrate quite similar behavior of the T W and their dependence 
upon the parameters (for more details, see [19]). 




6 Final remarks 

So it was shown in this study, that the generalized convection-reaction-diffusion equa- 
tion (JTJ) possesses a wide variety of TW solution, such as solitons, compactons, shock- 
like solutions and peakons. The result of the qualitative analysis and numerical simu- 
lation are in agreement with each other. The hyperbolic case occurs to be more reach, 
since the soliton-like solutions and peakons are observed only when a > 0. Surpris- 
ingly enough, one-sided compactons exist not only when the stationary point lies on 
the singular line A(U) = 0, but also when it lies left-hand side of it. Condition m > 1 
assures in the last case the existence of the sharp front. 

The results obtained are not rigorous. Only heuristic arguments are attached in 
favor of the statement that for m ^ 1/2 an infinite "time" is needed to reach the point 
(Uq, 0), moving along the incoming separatrice of the saddle sector. The precision of 
the numerical method used does not enable us to get convinced that we deal with shock 
fronts with exponentially localized "tails" in cases of small m G N + , when the critical 
point (Uq, 0) demonstrates essential instability and the "trains" of sharply ended pulses 
are observed, instead of truly localized solitary waves. More precise description of the 
wave patterns seems to be possible by presenting them in the form of exponential series 
|13] . Another possibility is connected with the application of the rigorous computing 
methods [23J. Aside of the scope of this study remained the study of stability and 
attracting properties of the invariant TW solutions. We plan to address these issues 
in the forthcoming studies. 
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